Geometric Entanglement in a One-Dimensional Valence Bond Solid State 
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In this paper we provide the analytical derivation of the global geometric entanglement per block 
for the valence bond solid ground state of the spin-1 AKLT chain. In particular, we show that 
this quantity saturates exponentially fast to a constant when the sizes of the blocks are sufficiently 
large. Our result provides the first known example of an analytical calculation of the geometric 
entanglement for a gapped quantum many-body system in one dimension and far away from a 
quantum critical point. 
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Introduction.- The study of entanglement in quantum 
many-body systems has attracted great attention in re- 
cent years. Quite notoriously, the establishment of a 
relation between entanglement and the Rcnormalization 
Group (RG) has proven particularly fruitful [1, 2, 3, 4]. 
In the case of a quantum system in one spatial dimension 
(ID) with a field theory limit, the c-thcorem of Zamolod- 
chikov [2] involved the existence of a global entanglement 
loss along successive RG transformations. And not long 
ago, it was proven that this loss of entanglement is ac- 
tually fine-grained, in the sense that it can be casted in 
a set of majorization relations for the reduced density 
matrices of the ground state of the system [3, 4]. 

Also for quantum systems in ID, it was recently shown 
how to quantify the idea of entanglement loss along RG 
flows in terms of the actual distance between the ground 
state of the system and the closest separable state in 
the Hilbert space. This novel idea was formulated in 
terms of the so-called global geometric entanglement per 
block, which was obtained for ID systems close to and 
at quantum critical points with an underlying conformal 
fleld theory [5, (j]. Despite of its apparent simplicity, 
the global geometric entanglement has proven very hard 
to compute. Regarding this, the results from Ref. [G] 
were significant since they managed to obtain analyti- 
cally this measure of entanglement for a whole class of 
extended quantum systems as long as they were close to 
criticality. Nevertheless, the studies of gapped systems 
in ID far away from a quantum critical point are almost 
non-existent and have always been restricted to numeri- 
cal calculations [7]. Exact and analytic examples of the 
behavior of the global geometric entanglement far away 
from ID quantum critical points are still missing. 

In this paper we address the above situation by com- 
puting analytically the global geometric entanglement per 
block for the ground state of the AKLT model. The AKLT 
model was introduced by Affleck, Kennedy, Lieb and 
Tasaki [8, 9], and constitutes the paradigm of a gapped 
spin-1 chain with Heisenberg-like interactions supporting 
Haldane's conjecture [10]. It has also the nice property 
that its ground state is a one-dimensional Valence Bond 
Solid (VBS) state which, in fact, is a particular case of a 



Matrix Product State (MPS) [11]. Our main result can 
be stated as follows: in the ground state of the ID spin- 
1 AKLT model, the global geometric entanglement per 
block of size L (which we call £{L)) obeys the law 



f(L) =log2-log 



1 



(1) 



in the thermodynamic limit, which saturates at large L 
to £{L > 1) - log 2. 

The AKLT model- The AKLT model was introduced 
by Affleck, Kennedy, Lieb and Tasaki in Refs. [8, 9] and 
has become a model of reference in condensed matter 
physics. Its importance relied originally in that it was 
the first analytical example of a quantum system that 
supported Haldane's conjecture [10], since in the ther- 
modynamic limit it is a local spin-1 Hamiltonian with 
Heisenberg interactions and a non-vanishing spin gap. 
More recently, it was shown how to simulate this model 
with cold atoms in optical lattices [13], and its utility for 
quantum computation was considered in [14]. Its ground 
state is a VBS state and, as such, it is closely related to 
the Laughlin wave function [ l ."i] and the fractional quan- 
tum HaU effect [ ; (i]. 

In the case of a quantum chain of N spins 1 in the bulk 
and two spins 1/2 at the boundaries, the Hamiltonian of 
the model reads 

r— 1 ^ ''^ 

(2) 

where S^^^ is a vector of spin-1 operators at site [r], and 
operators n^^'^^ and ttI^'^+^I are respectively given by 



,[0.1] 



1" 



[N.N+l] 



1'^ 



(3) 

where st°l and st^+^l are vectors of spin- 1/2 operators at 
sites [0] and [A^ -I- 1] . The ground state of the system can 
be represented in different ways. Here we choose to work 
with the representation 

!*> = |*-)["''ll|*-)[l'-2]...|^-)[A^',JV+l] ^ 

(4) 



2 



where = (|01) - |10))/V2 is a singlet state 

for a pair of ancillary spins 1/2 at sites [r',r + 1], and 
W^^'^' 1 is an isometry that projects a state of two spins 
1/2 on their symmetric subspace, which describes a state 
for a spin 1. The state constructed in this way is a 
one-dimensional VBS state, and can be proven to be 
the unique ground state of the Hamiltonian given in 
Eq. (2). Some of the entanglement properties of this 
state have been previously considered in the literature 
[16, 17, 18, 19, 20]. 

Geometric entanglement per block.- Let us now intro- 
duce the measure of entanglement that we shall use as 
follows. Consider a pure and possibly not normalized 
quantum state of TV parties l^*) <^H = "^''^ > where 

"HW is the Hilbert space of party i. As explained in Ref. 
[21], the global multipartite entanglement of j^*) can be 
quantified by considering the maximum fidelity [Amaxj 
between the quantum state and all the possible sep- 
arable and normalized states |<I>) of the N parties. 



(5) 



In order to have a well defined measure of entanglement 
we take the natural logarithm, 



^(*) = -log(|A, 



(6) 



We will be interested in the above quantity per party, 
which has a well defined limit when N ^ oo: 



En - N-^E{^) , 



E — lim En . 

N—*oo 



(7) 



The above quantity has been considered in a variety of 
situations [r), G, 7, 21, 22]. In this work we will always 
consider the case in which the one-dimcnsional quantum 
system is of infinite size. Moreover, the different parties 
will be contiguous blocks of L spins 1 as shown in Fig. 
(1). The entanglement measure E corresponds then to 
the global geometric entanglement per block of size L in 
the thermodynamic limit [o, G]. 



singlets in Eq. (4) as J^' ) = Y.a' ,ada' ,a\a' a) , where 
both indices a' and a can take the values or 1, and 
where do,i = — ^i^o = l/\/2- Also, the isometrics [''■'' 1 
in Eq. (4) are chosen such that they are equal at all 
sites ['"''' 1 = W V[r, r']. These can be written as 
W = Ej.a.a' ^a,a'l«)(aa'l' ^hcrc i = -1,0,1 and both 
a and a' can be either or 1. Importantly, in this work 
we choose W such that 



(8) 



where I3 is the identity operator in the three-dimensional 
Hilbert space of spin 1. This choice of W makes state \^) 
in Eq. (4) not to be normalized to one. 
Let us now write the VBS state as 



where indices Oq = 0, 1 and on+i =0,1 represent the 
states of the spins 1/2 at the boundaries, and indices 
V = — 1, 0, 1 for r = 1, • • • , iV represent the states of the 
N spins 1 in the bulk. According to the considerations 
above, coefficient Caj^.i^,... .i„.ajv+i can be decomposed as 

{ai-a'„} 

(10) 

Next, we decompose the coefficients da', a of the singlets in 
terms of the product of two matrices P and Q as da' .a = 
Pa Qa^ where index a can take the values or 1, and 
where the only non-zero components of P° and are 
given by 

r,l/4 ' ^ 1 91/4 



21/4 ' 1 21/4 



21/4 



(11) 



Considering together Eqs. (9), (10) and (11), our final 
MPS decomposition for the VBS state ]*) of Eq. (4) can 
be written in terms of matrices P and Q for the spins 
1/2 at the left and right boundaries respectively, and a 
tensor M of components 
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FIG. 1: (Color online) The ID spin chain is divided into con- 
tiguous blocks of L spins 1. For instance, in this diagram 
the system has been divided into contiguous blocks of L = 4 
spins. 

MPS representation of the ground state.- For the pur- 
pose of this paper, we represent the VBS state in Eq. 
(4) in terms of a MPS. Let us start by writing the 



(12) 



for the spins 1 in the bulk. This MPS is represented 
diagrammatically in Fig. (2.i). 

Geometric entanglement in the ground state.- As a first 
step towards the derivation of the global geometric en- 
tanglement per block, we compute the MPS transfer ma- 
trix A(\) of the ground state of the system. This trans- 
fer matrix, A{1) = Ea,/3,x5 ^(l)(7.f)V/3)(7^l (where 
(•] = l')"^), has components 



(13) 
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FIG. 2: (Color online) Diagrammatic representation of (i) an 
example of MPS decomposition of the VBS state in Eq. (4) 
for a system with three spins 1 in the bulk and two spins 1 /2 
at the boundaries, and (ii) MPS transfer matrix ^4(1) in Eq. 
(13). In these diagrams, tensors are represented by geometric 
shapes, and indices are represented by emergent legs. A leg 
shared by two shapes corresponds to the contraction of an 
index between two tensors, such that there is a sum over all 
its possible values. 



which are represented diagrammatically in Fig. (2.ii). 
Using Eqs. (8), (11), (12) and (13), we obtain the follow- 
ing block-diagonal expression for ^(1): 

A{1) =i(|ll)(00| + |00)(ll|) 
+ i(|ll)(ll| + |00)(00|-|10)(10|-|10)(01|). 



Diagonalizing the above matrix, we get its spectral de- 
composition 

^(1) = ^|ei)(e.i| - \ (|e2)(e2| + |e3)(e3| + |e4)(e4|) , 



(14) 

^(100) -111)), 



where |ei) = -i= (|00) + |11)), jea) 
leg) = |10) and \ei) - |01). 

Using the above decomposition we can now easily com- 
pute the norm (^'|^') of the state. For a system with open 
boundary conditions (OBC), where two spin-1/2 parti- 
cles are attached at the boundaries of the system (sites 
[0] and [TV -I- 1]), this norm is given by 



m^) = {L\A{ir\R)^ 



N 



(OBC) , (15) 



where the left and right vectors \L) and \R) are given 
by \L) ^ \R) = (|00) + 111)) /V2. Similarly, for periodic 
boundary conditions (PEC) (where all the particles have 



spin 1) the norm is given by 



(^l^-) ==tr {A{\ 



N 



N 



(PBC) 



(16) 

which matches the result for OBC in the limit N ^ oo, 
as expected. 

We are now in position to derive the global geometric 
entanglement per block of size L in the thermodynamic 
limit. This evaluation follows similar steps to the one in 
Ref. [G]. More precisely, our aim now is to compute the 
quantity 



|A| 



2_ (V1/|$)($|V1/) 



(17) 



maximized over all product states | $) of blocks of L con- 
tiguous spins 1. In particular, for i = 1 we have that 
I $) is a separable state of all the spins 1 that define the 
AKLT model in the thermodynamic limit. 

As explained in Ref. [(>], the maximization of Eq. (17) 
for a MPS reduces to the problem of finding 



Idl"^ = max\{r\* {r\A{L)\r)\ry 



(18) 



where A{L) is the MPS transfer matrix for blocks of 
length L, and the maximization is done over the vec- 
tor jr) with the normalization constraint (r|* • |r) = 1. 
Following Ref. [d], for ^ 1 we then have that 



> N 



lAI- 



f \d\' 



V(3/4)^ 



N 



(19) 



where L is the size of the block. 

In order to evaluate Eq. (18), we remember the prop- 
erty that A{L) = A{1)^ (see the discussion around Fig. 
(3) in Ref. [6] for a derivation of this property). Thanks 
to this fact, together with the spectral decomposition 
from Eq. (14) for A{1), we can evaluate A{L) for any 
block size L. The maximization in Eq. (18) can then be 
done in the following way: assume that |r) = a|0) -|-/3|1), 
where a and f3 are complex numbers. Then, we have that 



\r)\rr 



a 



^|00) + 1/3^111) + a*/3|01) 4- a/3*\10) . (20) 



In terms of a and (3 the quantity |c?p reads 

L / -, \ 



\d\' 



1 



a- 




(21) 



Since the vector |r) must be normalized to one, we have 



that 



1, and therefore 




(22) 



regardless of the values of a and /3. Quite conveniently, 
no maximization is needed at all in Eq. (18) [2G]. 
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Finally, using Eqs. (6), (7), (19) and (22) we get our 
result for the global geometric entanglement per block of 
size L in the thermodynamic limit, which reads 

£:(L)=log2-log|^l+ (^-i) j (23) 

as claimed in the introduction. 

As expected for a gapped quantum many-body system, 
Eq. (23) saturates for large block sizes L m£{L ^ oo) = 
log 2 ^ 0.693147. This saturation takes place exponen- 
tially fast and, as seen in Table (I), happens for relatively 
small values of L. Furthermore, note that £{L) decays to 



L 


£{L) 


L 


£(L) 


1 


1.098610 


8 


0.692995 


2 


0.587787 


9 


0.693198 


3 


0.730888 


10 


0.693130 


4 


0.680877 


11 


0.693153 


5 


0.697271 


12 


0.693145 


6 


0.691776 


13 


0.693148 


7 


0.693605 


14 


0.693147 



TABLE I: Some numerical values of £{L) for different block 
sizes L. At L = 14 we have that £{L = 14) ~ log 2 ~ 0.693147 
within the considered accuracy. 

a constant value equally fast as other entanglement mea- 
sures of the system, such as the entanglement entropy 
S{L) of a block of size L [17]: 

Notice also that both £{L) and S{L) decay as fast as 
the two-point correlation function for large L [9, IG], 

(S'^^'^S^^^) ^ (— l)^- A similar behavior is expected for 
the ground state of the more general bilinear biquadratic 
Hamihonian H = cos6'5'M^[''+il + sin6l(5M5['^+il)2 
(which reduces to Eq. (2) for 9 = 0.10247r) in the so- 
called Haldane phase -0.257r < 6 < 0.257r [2'>]. 

Conclusions.- Here we have computed analytically the 
global geometric entanglement per block for the valence 
bond solid ground state of the spin-1 AKLT chain in the 
thermodynamic limit. This quantity saturates to a con- 
stant exponentially fast for large block sizes. Our result 
is the first example of an analytical calculation of the ge- 
ometric entanglement for a gapped quantum many-body 
system in ID far away from a quantum critical point. 

The work in this paper opens the possibility of future 
analytical studies of the geometric entanglement in other 
gapped ID quantum many-body systems. For instance. 



it would be possible to apply the techniques of this pa- 
per to investigate the spin-S AKLT chain [19], the inho- 
mogeneous AKLT model [23] and the SU{n) generalized 
valence bond solid states [24]. 
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